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I .  INTRODUCTION 


A.  Obj ectives 

This  study  had  two  objectives.  The  first  was  to  assess  the  ability  of 
the  three-dimensional  (3-D)  HULL1*2  hydrocode  to  accurately  predict  shock 
diffraction  loading  on  all  faces  of  a  rectangular  parallelepiped  for  shock 
strengths  of  interest  in  the  study  of  the  airblast  vulnerability  of  tactical 
equipment.  A  second  objective  was  to  provide  further  quantification  of  the 
enhancement  in  reflected  overpressure  that  can  occur  for  a  shock  with  an 
overpressure  <  140  kPa  (20.3  psi)  when  it  strikes  a  target  at  an  angle  of 
obliquity3"6  (typically  near  45  degrees).  It  is  important  from  a  vulnerability 
viewpoint  to  determine  whether  or  not  such  an  enhanced  peak  can  have  a  suf¬ 
ficient  duration  and  spatial  extent  that  a  target  struck  at  obliquity  is  more 
vulnerable  than  one  struck  at  normal  incidence. 

A  combined  experimental  and  computational  research  program  was  undertaken 
at  the  Ballistic  Research  Laboratory  (BRL)  to  provide  these  assessments. 

This  report  documents  the  first  phase  of  the  program, 

B.  Background 

It  has  generally  been  assumed  that  the  most  severe  shock  loading  situa¬ 
tion  for  a  given  structure  occurs  when  the  shock  strikes  its  most  vulnerable 
face  at  normal  incidence.  However,  when  a  shock  wave  with  an  overpressure 
<  140  kPa  strikes  a  target  at  oblique  incidence,  the  peak  reflected  over¬ 
pressure  can  be  higher  than  that  for  normal  reflection. 3-6  This  oblique 


JJV.  A.  Fry,  ft.  E .  Durrett,  G.  P.  Ganong,  D .  d*  Matuska,  N.  D.  Stuoker, 

B,  S,  Chambers,  C.  E.  Needham,  and  C,  D.  Westmoreland,  "21 he  NULL 
Hydrodynamics  Computer  Code,"  AFVL-TR-76-183,  US  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Foroe  Base,  NM,  September  1976.  (AD  HB014070L) 

O 

“N.  A.  Fry,  C,  B.  Needham ,  ,V.  Stuoker,  B.  S.  Chamber  a.  III,  and 
G.  P.  Ganong,  "AFWL  NULL  Caloulatione  of  Air  Blast  Over  a  Dam  Slope ,  " 
AFWL-TR-76-154,  US  Air  Foroe  Weapons  Laboratory ,  Kirtland  Air  Foroe 
Base,  NH,  October  1976 .  (AD  §B016289L) 

3 

J.  von  Neumann,  " Oblique  Reflection  of  Shooks, "  Bureau  of  Ordnance 
Explosive  Research,  Report  12,  1943. 

Polaaheok  and  R.  J.  Seeger,  "Regular  Reflection  of  Shooks  in  Ideal 
Gases,"  Bureau  of  Ordnance  Explosives  Research,  Report  13,  1944. 

5 

L.  G.  Smith,  " Photographic  Investigation  of  the  Reflection  of  Plane  Shocks 
in  Air,"  Office  of  Scientific  Research  and  Development ,  Report  6271,  1945. 

j* 

'B.  P.  Bertrand,  "Measurements  of  Weak  Shock  Wave  Reflected  Pressure 
Nietoriee  on  a  ^^Dimensional  Surface , "  ABRL-MR-02966 ,  US  Amy  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  MD  October  1979.  (AD  §A 080S39) 


9 


interaction  process  is  of  interest  fToin  both  fluid  dynamic  and  vulnerability 
viewpoints.  Both  the  duration  and  magnitude  of  the  peak  reflected  overpres¬ 
sures  are  functions  of  the  incident  shock  overpressure  and  the  angle  of 
incidence  between  the  shock  wave  and  the  target  face.  The  duration  is  also  a 
function  of  the  distance  along  the  target  face,  measured  from  the  leading  edge. 
When  the  angle  of  incidence  is  o^,  the  angle  at  which  Mach  reflection  begins, 

the  reflected  overpressure  reaches  its  highest  possible  value.  For  relatively 
small  targets,  the  enhanced  peak  reflected  overpressure  is  difficult  to 
measure  experimentally  because  of  its  small  spatial  extent  and  short  time 
duration,  and  the  response  limitations  of  pressure  gages.  The  peak  reflected 
overpressure  is  difficult  to  predict  with  a  finite-difference  hydrodynamic 
computer  code  because  of  the  general  tendency  of  such  codes,  especially 
Eulerian  codes,  to  smear  discontinuities  such  as  the  incident  and  reflected 
shock  waves. 

The  peak  overpressure  enhancement  at  obliquity  for  shock  waves  with  an 
overpressure  <  140  kPa  is  of  interest  because  such  shock  waves  cover  a  large 
part  of  the  assumed  threat  range  from  tactical  nuclear  weapons.  The  S-280 
Electrical  Equipment  Shelter  (henceforth  called  the  shelter)  houses  commun¬ 
ications  equipment  for  many  Army  systems,  and  will  be  on  the  tactical  battle¬ 
field  in  large  numbers.  The  shelter  is  currently  being  hardened  to  decrease 
its  blast/thermal  vulnerability.  A  34. S  kPa  (S.O  psi)  overpressure  shock 
wave  was  chosen  for  the  present  work  because  it  represents  a  mid-range  threat 
level.  At  normal  incidence,  the  peak  reflected  overpressure  for  this  shock 
wave  is  78.5  kPa  (11.4  psi).  Mach  reflection  for  a  34,5  kPa  shock  wave  begins 
at  52.5  degrees  obliquity.  At  that  angle,  the  peak  reflected  overpressure  is 
estimated6  to  be  as  high  as  122  kPa  (17,7  psi),  55  percent  greater  than  the 
value  for  normal  incidence. 


C.  Topics  Discussed 

A  brief  introduction  to  the  principal  features  of  oblique  shock  diffrac¬ 
tion  is  provided,  with  particular  emphasis  on  the  circumstances  under  which 
the  peak  reflected  overpressure  can  exceed  that  for  normal  reflection.  The 
experiments  performed  are  described,  including  the  facility,  the  model, 
pressure  gage  types  and  locations,  and  the  shock  and  ambient  conditions.  The 
matching  hydrodynamic  computer  code  computations  are  also  described,  including 
detailed  descriptions  of  the  computational  grids  used. 

In  this  report,  the  computational  and  experimental  results  for  overpres¬ 
sure  are  compared  at  representative  points  on  each  face  of  the  shelter  model. 
Detailed  analyses  of  these  results  are  included  in  the  presentation  of  the 
data  for  each  face.  As  the  experiments  and  computations  progressed,  it 
became  apparent  that  neither  the  hydrocodc  computations  nor  the  experiments 
were  producing  computed  or  measured  overpressure  peaks  as  high  as  expected. 

The  study  was  then  broadened  to  include  c.  grid- resolution/convergence  analysis 
with  the  hydrocode,  and  a  gage-diameter/frequency-response  analysis  with  the 
experiments.  The  results  of  this  study  are  also  presented. 


II.  OBLIQUE  SHOCK  DIFFRACTION 


A.  Principal  Features 

When  an  incident  shock  wave  strikes  a  given  face  of  a  rectangular 
parallelepiped  at  some  oblique  angle,  it  also  strikes  another  face  at  the 
complement  of  that  angle.  The  initial  contact  between  the  incident  shock 
wave  and  the  structure  occurs  at  the  comer  formed  by  the  intersection  of 
these  two  faces.  This  comer  becomes  a  shock  diffraction  comer. 


As  the  divided  incident  shock  travels  along  the  two  perpendicular  faces, 
rarefaction  waves  emanating  from  the  diffraction  comer  travel  along  the  faces 
at  the  local  speed  of  sound  relative  to  the  moving  gas.  Depending  on  the 
angle  of  incidence  of  the  shock  on  the  face,  the  leading  edge  of  the  rare¬ 
faction  wave  will  either  fall  continually  farther  behind  the  incident  shock/ 
surface  intersection  point,  or  it  will  travel  coincident  with  it.  The  initial 
reflected  overpressure  experienced  at  a  given  point  on  the  surface  is  not 
relieved  until  the  comer  rarefaction  wave  arrives,  or  similar  waves  arrive 
from  other  parts  of  the  flow  field.  When  the  rarefaction  wave  is  traveling 
coincident  with  the  shock/surface  intersection  point,  the  reflected  overpressure 
is  relieved  immediately.  In  the  absence  of  viscous  effects,  this  results  in 
a  pressure  peak  of  vanishingly  small  spatial  extont  and  duration.  A  more 
detailed  discussion  of  oblique  shock  diffraction  may  be  found  in  Appendix  A. 


B.  Reflected  Shock  Enhancement 


The  incident  shock  wave  of  interest  here  is  a  34.5  kPa  overpressure  shock. 
Figure  1  shows  the  shock  reflection  factor  (the  ratio  of  the  peak  reflected 
overpressure  to  the  incident  overpressure)  as  a  function  of  the  angle  of  shock 
incidence,  a,  for  this  shock  strength.  Regular  reflection  theory  is  valid  for 
0  <  ct  <  «*£.  This  angle,  is  the  smallest  angle  at  which  the  comer  rare¬ 
faction  wave  travels  coincident  with  the  incident  shock/surfaee  intersection 
point.  The  experimental  data  indicate  that  the  peak  overpressure  is  reached 
at  w  S2.S  degrees.  At  that  angj .  a^,  the  intersection  point  between  the 

incident  and  reflected  shocks  is  on  the  verge  of  lifting  off  the  surface;  for 
o  >  a  Mach  stem  and  triple  point  are  formed. 


Information  inferred6  by  measuring  the  Mach  stem  velocity  for  a  ~ 

indicates  that  the  peak  reflected  overpressure  behind  the  Mach  stem  may  be 
even  larger  {122  kPa)  than  that  indicated  by  the  experimental  data  in  Figure  1. 
The  enhanced  peak  reflected  pressure  for  a  ~  is  of  interest  fro*  computational 

and  experimental  viewpoints,  and  may  have  important  implications  in  blast 
vulnerability. 


of  incidence.  (Reproduced  in  part  from  Reference  6) 


III.  EXPERIMENTS 


A.  Experimental  Facility  and  Model 

The  experiments  were  conducted  in  the  BRL  0.6  metre  diameter  shock  tube.7 
A  non-responding  1/18.45  scale  model  of  the-  shelter  was  constructed  of  aluminum, 
with  dimensions  19.58  cm  (width)  by  11.38  cm  (height)  by  11.79  cm  (depth). 

One  of  the  19.58  cm  by  11.38  cm  faces  is  defined  as  the  front  face.  The  model 
was  mounted  in  the  shock  tube  so  that  the  angle,  a,  between  the  front  face  and 
the  incident  shock  wave  front  was  52.5  ±  0.5  degrees.  One  of  the  11.79  cm  by 
11.38  cm  faces  is  defined  as  the  windward  side  face,  the  angle  between  that 
face  and  the  incident  shock  being  37.5  degrees.  The  model  produces  a  blockage 
of  9.4  percent  of  the  50.8  cm  by  50.8  cm  test  section  cross-sectional  area  in 
that ^orientation,  based  on  its  projected  area. 

B.  Experiments  Conducted 

Airblast  loading  experiments  were  conducted  for  various  overpressures  and 
angles  of  incidence.  The  two  described  in  Table  1  will  be  discussed  here. 

The  angle  between  the  shock  front  and  the  shelter  front  face  is  a,  p  is 
absolute  pressure,  T  is  temperature,  and  E.l  and  E.2  are  the  (simplified)  shot 
numbers  used  in  this  report.  Additional  information  on  the  experiments  is 
given  in  Appendix  B. 


TABLE  1.  EXPERIMENTAL  SHOTS 


Shot  Number 

a 

Ambient  Values 

Shock  Overpressure 

Actual  This  Report 

(Degrees) 

p(kPa) 

T(°C) 

(kPa) 

24-7^-126  E.l 

52.5 

101.42 

22.25 

32.75 

24-79-134  E.2 

52.5 

101.90 

24.31 

34.50 

C.  Pressure  Gages 

There  were  a  total  of  eight  pressure  gages  on  the  model,  mounted  flush  to 
the  model  surface.  Tablo  2  shows  the  gage  positions  where  pressure  measurements 
will  be  compared  with  the  hydrocode  computations.  The  gage  locations  are 
given  in  a  primed  coordinate  system,  having  its  origin  at  the  bottom  comer 
of  the  leading  vertical  edge  of  the  model.  The  front  and  back  faces  are 
constant  X*  planes,  the  side  faces  constant  Y',  and  the  top  and  bottom  faces 
constant  3*.  The  gage  positions  are  defined  by  a  letter-and-nuraber  pair;  the 
letter  denotes  the  face  and  the  number  denotes  the  gage  position  on  the  face. 


y 

G.,  A .  Coulter  and  B .  P.  Bertrand,  " BRL  Shook  Tube  Facility  for  the 
Simulation  of  Air  Blast  Effects, n  BRlr-MR-1685,  US  Amy  Ballietio 
Research  Laboratory,  Aberdeen  Proving  Ground,  UD,  August  196 5 . 

(AD  #475669) 


if 


All  gages  in  shots  E.l  and  E.2  had  a  sensitive  element  diameter  of  0.51  cm 
(0.20  in),  except  gage  F2  in  shot  E.2,  which  had  a  smaller  diameter  of  0.16  cm 
(0.063  in)  and  a  higher  frequency  response  than  the  other  gages.  (Reference 
8  incorrectly  described  the  diameter  of  this  latter  gage  as  0.32  cm.) 

TABLE  2.  GAGE  POSITIONS 


Gage  Number 

Face 

X’  (cm) 

Y*  (cm) 

Z'  (cm) 

F2* 

Front 

0.00 

14.68 

2.84 

F4 

Front 

0.00 

4.90 

8.53 

B2** 

Back 

11.79 

14.68 

2.84 

B4 

Back 

11.79 

4.90 

8.53 

SI 

Leeward  Side 

5.89 

19.58 

5.69 

S2 

Windward  Side 

5.89 

0.00 

5.69 

T1 

Top 

5.89 

14.68 

11.38 

T2 

Top 

5.89 

4.90 

11.38 

♦High  frequency  gage,  shot  E.2. 
♦♦No  gage,  shot  E.2. 


IV.  FINITE  DIFFERENCE  COMPUTATIONS 
A.  Hydrodynamic  Computer  Code 

The  airblast  version  of  the  HULL1*2  hydrodynamic  computer  code  was  used 
>  for  the  shock  diffraction  computations.  The  HULL  code  currently  in  use  at 

the  BRL  is  AFWL  version  8,  received  in  September  1978,  with  modifications  made 


i 
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by  the  BRL,  and  by  SA1 9  under  contract  to  the  BRL.  Some  of  the  modifications 
were  necessary  to  run  the  code  on  the  BRL’s  CDC  7600;  others  were  necessary 
to  convert  from  using  a  SAIL10  preprocessor  to  using  a  CDC-UPDATE/SAI-POST 
preprocessing  system.  Other  modifications  were  added  to  allow  the  input  of 
an  off-angle  step  shock  through  any  combination  of  the  left,  bottom,  and  back 
(aft)  boundaries  of  a  3-D  computational  grid  and  the  left  and  bottom  boundaries 
of  a  2-D  Cartesian  grid. 

The  HULL  hydrocode  uses  an  explicit  time  step,  predictor-corrector  method 
similar  to  a  Lax-Wendroff 1 1  scneme  to  solve  the  inviscid  Euler  equations.  Each 
computational  step  is  performed  in  two  phases,  a  Lagrangian  phase  where  flow 
field  cells  perform  work  on  one  another,  and  an  Eulerian  phase  where  material 
is  fluxed  across  cell  boundaries  using  a  donor  cell  method. 

B.  Confutations  Performed 


Table  3  lists  the  HULL  hydrocode  computations  described  in  this  report. 
Computation  H.l,  which  matches  experiment  2.1  in  shock  strength  and  ambient 
conditions,  will  be  discussed  in  detail  in  this  and  the  next  section.  The 
remaining  HULL  computations,  H.2  -  H.7,  match  experiment  E.2.  Results  from 
H.2  -  H.7  will  be  discussed  in  the  section  titled  "Convergence  Study."  The 
Courant-Friedrichs-Lewy  number  was  0.5  for  all  computations;  artificial  vis¬ 
cosity  was  not  used. 


TABLE  3.  HULL  HYDROCOPE  COMPUTATIONS 


Calculation  Number 

Dimensions 

Cell  Sizes 
Relative  to  H.l 

§rar 

Overpressure 

CkPa) 

Actual 

This  Report 

37.126 

H.l 

3-D 

1.0 

32.75 

37.0 

H.2 

3-D 

1.0 

34.50 

37.001 

H.3 

3-D 

0.5 

34.50 

37.11 

H.4 

2-D 

1.0 

34.50 

37.9 

H.5 

2-D 

o.s 

34.50 

37.10 

H.6 

2-D 

0.25 

34.50 

37.6 

H.7 

2-D 

0.125 

34.50 

9J.  .4.  Hasdai,  8,  S.  Chatxtoers,  and  R.  W.  Clemens ,  "Suoport  to  BRL;  HULL  Code 
Implementation  on  a  CDC  *600, w  SAZ-8Q-701-AQ,  Science  Applications,  Inc,, 

NoLean,  VA,  August  1970* 

C.  Graham,  L.  P.  Gaby,  end  C,  E,  Rhodes,  C SAIL ,  A  x  Automated  Approach  to 
Software  Development  and  Management,  **  AWL  Interim  Report  1971-6,  US  Air 
Force  Weapons  Laboratory,  Klrtland  Air  Force  Base,  UN,  October  1976* 

*3i?.  D.  Riohtmyer  and  ff.  W.  Norton,  Difference  Methods  for  Initial  Value  Problem, 
Intersoisnce  Publishers,  Tea.,  John  Wiley  &  Sons,  Inc Second  Edition,  1967* 
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C.  Finite  Difference  Grids 


The  finite  difference  grid  used  for  computation  H.l  contains  92,512 
flow  field  cells,  with  a  49  *  59  *  32  grid  in  the  X,  Y,  and  Z  directions, 
respectivoly.  The  shelter  is  modeled  using  6,912  nearly  cubical  rigid  cells, 
with  16  equal  X-direction  cells  (AX  =  0.7366  cm),  27  equal  Y-direction  cells 
(AY  =  0.7253  cm),  and  16  equal  Z-direction  cells  (AZ  =  0.7112  cm).  The 
rigid  cells  do  not  directly  enter  the  computation,  but  they  do  require  the 
same  storage  space  as  a  hydrodynamic  cell,  so  this  space  is  wasted.  To 
minimize  the  smearing  of  the  computational  shock  as  it  passes  through  the 
grid  prior  to  stril.ing  the  shelter,  the  shock  is  placed  well  into  the  grid 
at  the  initiation  of  the  computation,  2.66  cm  upstream  from  the  shelter  leading 
edge.  The  shock  input  algorithm  keeps  track  of  the  theoretical  intersection 
of  the  shock  wave  with  the  computational  boundaries,  progressively  moving 
the  input  shock  alo-  j  the  boundaries  as  the  computation  proceeds. 

Appendix  C  contains  detailed  descriptions  of  the  computational  grids 
used  in  this  study. 

D.  The  Flow  Field 


Figure  2  shows  a  top  view  of  an  isobar  (constant  pressure)  plot  of  the 
flow  field  for  H.l  in  the  bottom- most  plane  of  cells  after  one  computational 
cycle.  (A  similar  plot  is  not  ava-‘lable  at  cycle  0.)  Shock  arrival  time 
at  the  shelter  leading  edge  is  defined  as  t  =  0.0.  In  Figure  2,  the  incident 
computational  shock  is  indicated  by  the  closely-spaced  pressure  contours. 

(The  apparent,  excessive  width  of  the  shock  is  an  artifact  of  the  contour¬ 
plotting  algorithm.)  The  shcck  is  moving  from  the  lower  left  corner  of  the 
figure  toward  the  upper  right  corner,  hnat  ippears  to  be  a  bent  right  end  of 
the  shock  is  actually  the  set  of  pressure  contour  indentification  numbers 
which  are  overwritten  by  the  plot  routine.  The  contour  labelled  "1"  Is  not  a 
pressure  contour,  but  i?  an  artificial  use  of  the  contour  algorithm  to  show 
the  outline  of  the  shelter  in  this  piano. 

Figure  3  shows  isobars  in  the  same  j*’°ne  at  t  »  374.9  us;  the  shock 
wave  has  passed  slightly  more  than  half-way  across  the  shelter.  Contour  5 
shows  the  general  shape  of  the  reflected  shc:k.  Figures  2  and  3  give  a 
qualitative  indication  of  the  flow  fiold  predicted  by  the  3-D  HULL  hydrocodc . 
The  next  section  provides  uetailed  comparisons  betweon  experiment  E.l  and 
HULL  computation  H.l. 


V.  COMPARISON  OF  COMPUTATIONAL  AND  EXPERIMENTAL  RESULTS 
A.  General  Comments 


This  section  shows  typical  comparisons  between  measured  and  computed 
overpressure  on  each  face  of  the  model.  In  this  section,  ,,cxpcriment,'  refers 
to  oxpetiaent  E.l,  and  "computation'*  refers  to  computation  H.l.  Computation 
H.l  was  run  on  the  BRL  CDC  7600  for  155  computational  cycles,  with  a  total 
job  time  of  49  minutes. 
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Figure  3.  Isobars  at  t  *  374.9  ys  in  ground  piano  for  H.l. 


B.  Front  Face 


Figure  4  shows  a  comparison  between  the  measured  and  the  computed  over¬ 
pressure  vs  time  for  gage  F2.  The  gage  is  located  at  3/4  of  the  distance  along 
the  shelter  front  face  from  the  leading  edge  and  at  1/4  of  the  shelter  height 
from  the  ground  plane.  The  agreement  between  the  two  results  is  good,  except 
for  the  initial  shock  interaction  with  the  structure.  Because  the  computational 
shock  is  spread  over  three  or  four  cells,  it  shows  an  earlier  initial  rise  and  a 
reduced  peak.12  The  agreement  beyond  the  peak  is  very  good,  although  the 
computed  overpressure  is  slightly  greater  than  that  for  the  experiment  for 
0.6  <  t  <  1.2  ms.  The  second  peak  in  the  experimental  data  at  1.22  ms  is 
caused  by  the  arrival  of  a  secondary  shock  produced  by  the  interaction  of  the 
incident  shock  with  the  model;  it  traveled  to  the  shock  tube  wall,  reflected 
from  the  wall,  and  returned  to  strike  the  model.  A  similar  second  peak  is 
seen  in  the  other  comparisons  in  this  section.  The  X's  on  the  H.l  curve  in 
this  and  the  next  four  figures  mark  every  fifth  computed  data  point;  they  have 
been  added  primarily  as  a  visual  aid. 

Figure  4  shows  a  measured  peak  reflected  overpressure  of  86.2  kPa  occurring 
at  0.320  ms,  and  a  corresponding  computed  peak  of  76.2  kPa  at  0.331  ms.  The 
difference  in  time  is  due  in  part  to  the  difficulty  in  establishing  a  zero 
reference  time  at  the  leading  edge  of  the  shelter  for  the  experiment  because 
of  the  discrete  data  sampling  rate  and  the  gage  diameter  and  response  time. 

Time  zero  was  deduced  from  the  experimental  data  using  the  incident  shock 
speed  and  the  time  at  which  a  reference  gage  on  the  model  first  sensed  a  pres¬ 
sure  signal  above  a  threshold  value.  This  was  taken  to  indicate  shock  arrival 
at  the  gage,  but  this  method  has  some  inherent  uncertainties.  The  experimental 
wave  speed  was  computed  by  using  the  measured  pressure  jump  across  the  incident 
shock  wave.  The  wave  speed  thus  computed  for  B.l  was  389,3  m/s.  The  data  for 
E.l  were  sampled  at  2.5  us  intervals.  The  diameter  of  the  gage  was  0.51  cm, 
so  it  took  the  incident  shock  13.1  vs  or  5.24  sampling  counts  to  traverse  the 
gage.  The  time  difference  in  the  peaks  for  E.l  and  H.l  for  F.2  was  11  us, 
which  is  within  one  digital  sampling  count  of  the  13.1  us  shock  traversal 
time  across  the  gage. 

There  was  also  an  inherent  error  in  the  time  of  peak  overpressure  computed 
by  the  HULL  hydrocode.  The  numerical  diffusion  of  the  shock  wave  as  it  travels 
through  the  finite  difference  grid  has  been  demonstrated  12-14  to  cause  a 


IP  ♦ 

R.  £’.  Lottero,  "Comparison  of  3-D  Hydrocode  Computations  for  Shook  Diffraction 

Loading  on  an  8-280  Electrical  Equipment  Shelter,"  Army  Research  Office  Report 

80-3,  Proceedings  of  the  I860  Amy  numerical  Analysis  and  Computers  Conference, 

August  1980.  (AD  8A089089) 

*^R.  E,  Lottero,  "A  Detailed  Comparison  of  3-D  Hydro code  Computations  for  Shook 
Diffraction  Loading  of  an  S*  280  Electrical  Equipment  Shelter,"  ARBRL-TR-  02334,  US 
Amy  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground,  ND,  June  1981. 
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14 r,  a.  Gentry ,  L.  R.  Stein,  and  C.  V.  llirt,  "Three-Dimensional  Computer  Awlyois 
of  Shook  Loads  on  a  Simple  Structure,"  8RL-CR-219,  US  Amy  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  KD,  March  1978.  (AD  HB003208L) 
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Figure  4.  Comparison  of  measured  and  predicted  overpressure 
on  front  face  gage  P2. 


delayed  and  reduced  peak  reflected  overpressure  when  the  shock  interacts  with 
a  structure.  Thus,  each  method  contributed  to  the  uncertainty  in  establishing 
a  common  reference  time.  It  was  decided  to  establish  time-zero  estimates  for 
the  experiment  and  the  computation  independently  and  use  these  to  compare 
the  results,  rather  than  arbitrarily  shift  the  overpressure  curves  in  time  to 
get  a  more  pleasing  visual  alignment. 

Neither  the  measured  nor  the  computed  peak  reflected  overpressure  were 
near  the  expected  peak  of  115  kPa  (see  the  section  entitled  "Convergence 
Study").  The  measured  peak  of  86.2  kPa  is  25  percent  less  than  the  expected 
peak,  and  16  percent  greater  than  the  peak  for  normal  reflection.  The  computed 
peak  is  34  percent  less  than  the  expected  peak,  and  3  percent  greater  than 
the  peak  for  normal  reflection^  This  is  most  likely  due  to  the  finite  gage 
and  cell  sizes  which  are  used.  In  the  experiment,  the  gage  diameter  is  2.6 
percent  of  the  span  of  the  front  face.  In  the  computation,  the  cell  dimension 
in  the  Y  direction  on  the  front  face  is  42  percent  larger  than  the  gage  diameter 
and  3.7  percent  of  the  span.  It  appears  that  the  spatial  extent  of  the  over¬ 
pressure  peak  is  so  small  that  it  is  largely  integrated  out,  even  with  these 
relatively  small  units  of  measure.  Thus,  it  can  be  concluded  that  the  over¬ 
pressure  peak  does  not  make  a  significant  contribution  to  the  loading  on  the 
model  and  hence  is  probably  not  significant  to  the  shelter  itself.  The  problem 
of  resolving  the  peak  is  still  of  interest  from  a  fluid  dynamics  viewpoint,  and 
will  be  discussed  further  in  the  section  entitled  "Convergence  Study." 

Analysis  of  the  experimental  data  shows  that  gage  F2  recorded  a  higher 
peak  reflected  overpressure  (86.2  kPa)  than  did  gage  F4  (73.8  kPa).  This  is 
because  gage  F2  was  farther  from  the  leading  edge  than  was  gage  F4,  and  so 
the  pressure  peak  had  more  opportunity  to  develop  spatially,  as  discussed  in 
Appendix  A.  (Data  for  gage  position  F4  are  not  shown  here.) 

C.  Windward  Side  Face 


Figure  5  shows  a  comparison  between  the  measured  and  computed  overpressure 
histories  for  gage  S2,  located  in  the  center  of  the  windward  side  face.  The 
angle  between  the  shock  front  and  this  face  is  37.5  degrees.  The  agreement  is 
good  except  for  the  initial  shock  interaction.  For  this  angle,  the  expected 
peak  reflected  overpressure  (see  Figure  1)  is  essentially  equal  to  the  normal 
reflection  overpressure  of  74.2  kPa.  The  computational  peak  is  77,8  kPa,  5 
percent  greater  than  the  normal  reflection  value  and  7  percent  greater  than 
the  measured  peak  (73.0  kPa),  which  is  in  turn  2  percent  lower  than  the  normal 
reflection  peak.  The  peak  values  for  this  gage  position  will  be  discussed 
further  in  the  section  entitled  "Convergence  Study." 

D.  Leeward  Side  Face 


Figure  6  shows  a  similar  comparison  for  gage  SI,  located  at  the  center  of 
the  leeward  side  face.  It,  too,  shows  a  smearing  of  the  computational  shock, 
which  has  also  been  weakened  by  a  rarefaction  wave  produced  at  the  trailing 
edge  of  the  front  face.  The  general  agreement  between  the  curves  is  good. 

They  both  show  a  pressure  plateau  of  22.5  kPa  for  0.6  <  t  <  0.9  ms,  caused  by 
the  weakened  incident  shock.  The  pressure  rise  which  begins  at  0.9  ms  is 
caused  by  the  incident  shock  breaking  over  the  top  face,  sending  another 
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Figure  6.  Comparison  of  measured  and  predicted  overpressure 
on  leeward  side  face  gage  SI. 


weakened  shock  down  the  leeward  side  face.  The  curves  agree  qualitatively, 
but  the  computed  results  vary  from  5  to  10  percent  less  than  the  experimental 
results  after  0.9  ms.  It  may  be  that  viscous  effects,  which  HULL  cannot  model, 
are  becoming  important  by  this  time. 

E.  Back  Face 


Figure  7  shows  a  similar  comparison  for  back  face  gage  B2.  This  gage  is 
located  at  3/4  of  the  distance  along  the  back  face,  and  1/4  of  the  height  of 
the  shelter  from  the  ground  plane.  The  agreement  between  the  computed  and 
the  measured  results  is  very  good,  except  at  the  time  of  initial  shock  arrival 
and  for  t  >  1.0  ms. 

F.  Top  Face 

Figure  8  shows  a  comparison  between  the  measured  and  computed  overpressure 
histories  for  gage  T2,  located  on  the  top  face  (see  Table  2) .  Except  for  the 
overpressure  peaks,  the  agreement  between  the  computation  and  experiment  is 
good,  with  the  computed  pressure  consistently  larger.  The  computed  overpressure 
peak  of  39.5  kPa  at  0.22  ms  is  18  percent  greater  than  the  experimental  peak 
of  33.4  kPa  at  0.20  ms.  The  HULL  result  for  the  overpressure  peak  is  the 
value  in  error;  this  is  possibly  related  to  a  numerical  overshoot  associated 
with  the  top  corner  of  the  shelter.  The  comparison  plot  for  top  face  gage  T1 
(not  shown  here)  shows  similar  behavior. 


VI.  CONVERGENCE  STUDY 


A.  Reason  for  Study 

Because  the  expected  peak  reflected  overpressure  was  not  obtained  on  the 
front  face  (gage  position  F2)  in  either  tho  calculations  or  the  experiment, 
a  convergence  study  was  performed  to  study  the  sensitivity  of  the  peak  values 
to  both  cell  and  gage  size.  It  was  expected  that  smaller  cell  sizes  and  smaller 
gages  having  higher  frequency  response  were  needed  to  resolve  the  small  region 
of  enhanced  reflected  overpressure  on  the  front  faco. 

B.  Computations  and  Experiments 

Only  a  limited  grid  refinement  was  possible  for  the  3-D  computational 
problem  because  of  cost  and  storage  limitations.  Host  of  the  computational 
results  were  obtained  for  a  representative  2-D  slice  of  the  grid  for  11.1,  which 
then  modeled  a  shelter  with  infinite  height.  Experiment  E.2  was  performed  to 
study  the  effect  of  reduced  gage  size  and  higher  frequency  response  on  the 
measured  front  face  peak. 

The  cell-size  convergence  study  matching  shot  B.2  consisted  of  two  3-D 
computations,  H.2  and  H.3,  and  four  2-D  computations,  H.4  -  H.7  (sec  Table  3). 
The  relative  cell  sizes  in  Table  3  were  computed  by  dividing  the  cell  sizes 
for  the  particular  grid  by  the  respective  values  of  AX,  AY,  and  A2  on  the 
shelter  surface  for  H.l.  The  experimental  gage  sizes  are  scaled  in  the  same 
way,  using  AY  from  H.l  for  the  front  face  gage  and  AX  from  H.l  for  the 
windward  side  gage. 


24 


EXPERIMENT 

COMPUTATION 


C.  Results 


1.  Front  Face. 

Figure  9  shows  a  comparison  of  the  peak  reflected  overpressures?  obtained 
at  gage  position  F2  on  the  front  face.  Neither  the  experimental  (E.l)  nor  the 
computed  (H.l)  values  were  near  the  expected  peak  of  115  kPa  for  that  32.75  kPa 
shock.*  The  experimentally  measured  peak  of  86.2  kPa  is  25  percent  below  the 
expected  peak  and  16  percent  above  the  normal  reflection  peak  of  74.2  kPa. 

The  corresponding  computed  peak  in  H.l  is  76.2  kPa,  34  percent  less  than  the 
expected  peak  and  3  percent  greater  than  the  normal  reflection  peak. 

Once  these  comparisons  were  made,  the  problem  then  was  to  determine  what 
could  be  done  both  experimentally  and  computationally  to  obtain  peak  overpres¬ 
sures  closer  to  the  expected  values.  If  neither  E.l  nor  H.l  produced  peak 
overpressures  ,.<’ur  115  kPa  because  the  measuring  elements  (either  gages  or 
flow  field  cells)  were  too  large,  then  it  should  be  possible  to  measure  higher 
peaks  expe.  mentally  with  a  smaller  diameter,  higher-frequency-response  gage 
and  compute  higher  peaks  with  smaller  flow  field  cells.  For  E.l,  the  gage 
diameter  was  0.51  cm,  2.6  percent  of  the  span  of  the  front  face.  For  H.l,  the 
cell  dimension  in  the  Y  direction  on  the  front  face  was  0.7253  cm,  3.7  percent 
of  the  span  of  the  front  face,  and  42  percent  larger  than  the  gage  diameter. 

A  second  experiment,  shot  E.2,  was  fired  with  a  smaller  diameter  (0.16  cm) 
'.igher  frequency-response  gage  at  position  F2.  The  experimental  conditions 
and  gages  used  for  E.2  are  summarized  in  Tables  1  and  2.  The  measured  peak 
overpressure  at  gage  F2  in  shot  E.2  was  88.6  kPa,  27  percent  below  the  expected 
peak  of  122  kPa,  and  13  percent  above  the  normal  reflection  peak  of  78.5  kPa. 

The  relative  difference  in  measured  reaks  for  E.l  and  E.2  was  within  the  range 
of  experimental  error,  and  hence  inconclusive. 

The  diameter  of  the  gage  at  position  F2  for  shot  E.2  was  0.82  percent  of 
the  length  of  the  front  face,  a  reduction  in  gage  diameter  of  69  percent  relative 
to  shot  E.l.  Even  with  this  relatively  small  high-frequency-response  gage,  the 
overpressure  peak  was  of  sufficiently  short  duration  and  limited  spatial  extent 
that  it  was  of  little  significance  to  the  loading  on  the  shelter  front  face. 

The  computational  convergence  study  involved  several  computations,  both 
2-D  and  3-D,  for  a  34.5  kPa  overpressure  shock  (see  Table  3).  The  column  in 
Table  3  indicating  relative  cell  si«es  is  computed  by  taking  a  cell  dimension 
on  the  model  surface,  and  dividing  it  by  the  corresponding  cell  dimension  in 
computation  H.l.  (See  appendix  C  for  additional  details  concerning  the  compu¬ 
tational  grids.)  Table  4  gives  a  summary  of  the  computed  and  measured  peak 
overpressures  at  gage  position  F2.  (These  values  are  also  plotted  on  Figure  9.) 
There  is  a  clear  trend  toward  increasing  peak  overpressure  with  decreasing  flow  - 


*  The  expected  peak  reflected  ovorpressure  for  a  34.5  kPa  shock  at  52.5  degrees 
obliquity  is  122  kPa,6  1.55  times  the  normal  reflection  overpressure  of  78.5 
kPa.  If  a  similar  enhancement  is  assumed  for  the  32.75  kPa  overpressure  shock 
in  E.l,  then  a  peak  reflected  overpressure  no  higher  than  115  kPa  can  be 
expected.  This  is  because  aM  >  52.5  degrees  for  a  32.75  kPa  shock. 
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field  cell  size.  However,  because  of  the  rapidly  increasing  cost  of  the 
computations  as  the  cell  size  was  decreased,  it  was  not  possible  to  determine 
whether  the  computed  peak  would  approach  the  expected  value  asymptotically  as 
cell  size  was  decreased  further,  approach  some  other  value,  or  increase  with¬ 
out  bound. 


TABLE  4.  PEAK  OVERPRESSURE,  GAGE  F2 


Source 

Peak  Overpressure 
(kPa) 

Peak/NRP* 

Peak/IP** 

« 

86.2 

1.16 

.75 

H.l 

76.2 

1.50 

.66 

E.2 

88.6 

1.13 

.73 

H.2 

80.4 

1.02 

.66 

H.3 

89.4 

1.14 

.73 

H.4 

76.5 

.97 

.63 

H.S 

83.9 

1.07 

.69 

H.6 

92.2 

1.17 

.76 

H.  7 

100.4 

1.28 

.82 

*NRP  a  Theoretical  normal  reflection  peak  overpressure,  74.2  kPa  for  the 
32.75  kPa  overpressure  stock,  and  78.5  kPa  for  the  34,5  kPa  over¬ 
pressure  stock. 

**IP  «  inferred  peak  at  o^,  115  kPa  for  the  32.75  kPa  overpressure  stock, 
and  122  kPa  for  the  34.5  kPa  overpressure  shock. 

There  was  one  other  interesting  feature  noted  in  this  .convergence  study, 
iVft  directly  related  to  the  problem  at  hand.  There  is  a  systematic  displace¬ 
ment  between  the  2-D  and  3-D  results  for  identical  shock  and  ambient  conditions, 
apv*x  sly  due  to  <  inferences  in  the  finite  difference  algorithms  in  the  two 
version:-  of  the  HULL  code. 

2,  Windward  Side  Pace. 

Figure  10  shows  the  results  of  the  convergence  study  for  gage  position 
S2.  In  this  case,  a  «  37.5  degrees,  and  the  expected  peak  overpressures 
are  the  same  as  for  normal  reflection.  The  qualitative  trend  is  similar  to 


4  E.1 


Figure  10. 


Peek  overpressure  at  windward  side  face  gage  positio 


that  shown  in  Figure  9.  The  computational  results  overshoot  the  expected 
values  by  as  much  as  20  percent  as  the  cell  size  is  decreased,  indicating  an 
artificial  sensitivity  of  the-  peak  value  to  cell  size. 


In  this  study,  the  computed  peak  rises  monotonically  with  decreasing 
cell  size  for  both  the  front  and  windward-side  faces.  This  sensitivity  to 
grid  size  has  also  been  documented  by  Carpenter  et  al.15  for  strong  shocks 
(pressure  ratio  >  20),  but  the  relation  there  was  nonmonotonic. 


VII.  CONCLUSION 

This  study  has  shown  that  the  3-D  HULL  hydrocode  can  produce  accurate 
shock  diffraction  loading  predictions  for  this  class  of  problems  at  moderate 
cost.  The  BRL  code  version  has  been  modified  to  treat  an  incident  step  shock 
moving  obliquely  through  the  2-D  and  3-D  Cartesian  grids.  This  allows  the 
modeling  of  the  32.7,5  and  34.5  kPa  (nominal  5  psi)  overpressure  shocks  striking 
the  front  face  of  an  S-280  Electrical  Shelter  at  a  52.5  degree  angle  of 
incidence,  using  smooth  shelter  walls  instead  of  undesirable  rough  "stair- 
stepped”  surfaces  formed  if  the  shelter  were  rotated  in  the  grid.  The  agree¬ 
ment  between  the  3-D  computational  and  experimental  results  is  good  on  all 
shelter  faces.  Errors  of  from  -12  percent  to  -9  percent  (H.l/E.l  and  H.2/E.2), 
respectively)  are  experienced  in  resolving  the  peak  reflected  overpressure  on 
the  front  face  (P2),  and  +7  percent  to  +8  percent  on  the  windward  side  face 
(S2);  error  magnitudes  less  than  10  percent  are  present  beyond  the  peak.  A 
grid  convergence  study  has  quantified  the  sensitivity  of  the  peak  reflected 
overpressure  to  grid  size.  The  expected  enhanced  peak  overpressure  on  the 
front  face  at  52.5  degrees  obliquity  proved  to  be  difficult  to  measure  or 
compute.  It  is  of  such  short  duration  and  limited  spatial  extent  that  it  is 
probably  unimportant  as  an  enhanced-damage  mechanism  for  the  S-280  Electrical 
Equipment  Shelter. 
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APPENDIX  A 

OBLIQUE  SHOCK  DIFFRACTION 
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Figure  A-l  shows  a  schematic  view  of  oblique  shock  diffraction  over  a 
three-dimensional  box- like  target.  Figure  A- la  shows  the  target  sitting  on 
the  ground  with  the  line  labelled  "TARGET  AXIS"  normal  to  one  of  the  faces. 

An  incident  shock  is  propagating  toward  the  target,  with  the  shock  wave  front 
perpendicular  to  the  line  labelled  "SHOCK  PROPAGATION  AXIS,"  which  is  in  turn 
at  some  general  angle  to  the  target  axis.  Figure  A- lb  shows  the  shock  at  a 
later  time  when  the  shock  has  passed  part  of  the  way  across  the  target.  The 
two  vertical  faces  of  the  target  upon  which  the  incident  shock  impinges  may 
experience  different  pressure  loading  histories,  which  in  general  depend  on 
the  shock  strength  and  its  angle  of  incidence. 

Figure  A-2  shows  qualitatively  the  several  possible  wave  configurations. 
The  incident  shock  wave  front  is  denoted  by  S,  the  reflected  shock  wave  by 
SR,  the  rarefaction  wave  by  R,  the  wedge  corner  by  C,  the  intersection  point 

of  either  the  incident  shock  or  Mach  stem  with  the  wedge  by  T,  and  the  inter¬ 
section  of  the  rarefaction  wave  with  the  wedge  by  A.  The  side  of  the  wedge 
initially  facing  toward  the  shock  is  termed  side  1,  and  the  other  side  making 
a  right  angle  at  C  is  side  2.  Figure  A-2a  shows  the  wave  configuration  when 
the  incident  shock  wave  strikes  the  wedge  at  normal  incidence.  The  gas  behind 

S  is  brought  to  rest  on  side  1,  producing  reflected  shock  S_.  A  rarefaction 

R 

wave  R  is  generated  at  corner  C,  relieving  the  reflected  overpressure  behind 
SR.  On  side  2,  the  incident  shock  S  travels  past  the  corner  and  the  loading 

is  initially  unaffected  by  the  process  on  side  1  until  waves  generated  by  R 
move  around  corner  C  and  catch  up  to  S.  Figure  A- 2b  shows  the  wave  configur¬ 
ation  when  the  wedge  is  tilted  so  that  side  1  is  at  a  slight  angle  to  S.  The 
intersection  point  T  moves  away  from  comer  C  at  a  higher  speed  than  does 
point  A,  producing  three  distinct  regions  on  side  1.  The  region  beyond  T 
and  in  front  of  S  is  ambient  gas;  the  region  between  A  and  T  has  experienced 
only  the  oblique  reflection  of  S;  and  the  region  between  C  and  A  has  exper¬ 
ienced  not  only  the  oblique  reflection  of  S,  but  also  the  passage  of  rare¬ 
faction  wave  R.  The  action  of  S  on  side  2  is  that  of  a  weak  oblique 
reflection. 

Figure  A-2c  is  conceptually  the  same  as  Figure  A-2b.  The  angle  of  inci¬ 
dence  between  S  and  side  1  of  the  wedge  is  larger  in  Figure  A-2c  than  in 
Figure  A- 2b,  so  the  speed  of  T  along  side  1  had  decreased  and  the  distance 
between  A  and  T  now  is  growing  at  a  slower  rate  than  in  Figure  A-2b. 

Figure  A-2d  shows  the  situation  where  the  angle  between  S  and  side  1  has 
reached  the  catch-up  angle  a,,  where  the  rarefaction  wave  R  is  traveling  at 

exactly  the  same  speed  as  T.  In  this  case,  the  overpressure  due  to  the  shock 
reflection  is  relieved  immediately  for  these  weak  shocks.  Figure  A-2e  shows 
the  wave  configuration  for  an  angle  a  few  degrees  larger  than  ac;  the  inflected 

shock  SR  is  merging  with  the  incident  shock  S,  forming  a  Mach  stem  and  a  triple 

point  at  the  intersection  of  the  three  shocks.  The  angle  at  which  the  Mach 
stem  begins  to  fom  is  denoted  by  it  is  at  this  angle  that  the  peak 

reflected  overpressure  occurs.* 


•For  strong  shocks,  irregular  Mach  reflection  can  occur  where  a  compression 
wave  or  shock  can  follow  the  first  Mach  stem . 


(a)  Normal  reflection. 


(b)  Regular  reflection, 
near  normal. 


(c)  Regular  reflection 
near  <*c. 


Figure  A- 2.  Reflected  shock  configurations.  (Reproduced  from 
Reference  6.) 


Figure  1  (in  the  main  body  of  the  report)  shows  the  shock  reflection 
factor  (the  ratio  of  the  peak  reflected  overpressure  to  the  incident  over¬ 
pressure)  as  a  function  of  the  angle  of  shock  incidence,  ot,  for  a  34.5  kPa 
(5.0  psi)  overpressure  incident  shock  wave.  For  this  wave,  which  is  of  inter¬ 
est  in  this  study,  =  52.5°.  The  experimental  data  in  Figure  1  were  taken 

on  a  long  two-dimensional  wedge  by  Bertrand.6  As  may  be  seen  in  Figure  1, 
there  is  little  change  in  the  reflected  overpressure  from  a  =  0  (normal 
reflection)  to  a  =  36°.  Figure  1  does  show  that  for  values  near  o  =  52.5°, 
the  peak  overpressure  for  a  34.5  kPa  shock  is  considerably  larger  than  that 
for  normal  reflection.  This  peak  is  difficult  to  measure  experimentally  and 
difficult  to  predict  computationally  because  of  its  short  duration  and,  at 
least  on  small  models,  its  small  spatial  extent.  However,  for  large  buildings 
such  as  manufacturing  plants,  an  enhanced  peak  capable  of  causing  damage  may 
develop,  hence  the  interest 'in  oblique  reflection  of  low-overpressure  shocks. 


APPENDIX  B 

ADDITIONAL  INFORMATION  ON  EXPERIMENTS 


Table  B-l  shows  the  locations  of  the  centers  of  all  the  pressure  gage 
positions  on  the  S-280  shelter  model.  They  are  given  in  a  primed  coordinate 
system  relative  to  the  bottom  comer  of  the  leading  vertical  edge  of  the  model. 
The  leading  edge  is  defined  a.  the  vertical  edge  which  the  incident  shock  wave 
contacts  first.  The  bottom  comer  of  the  leading  edge  is  defined  as  (X*  =  0.0, 
Y’  =  0.0,  V  =  0.0),  where  X'  =  depth,  Y'  =  width,  and  Z'  =  height.  The  primed 
coordinate  system  is  defined  so  that  the  front  and  back  faces  of  the  model  are 
planes  of  constant  Xr;  the  side  faces  are  planes  of  constant  Y  * ;  and  the  top 
and  bottom  faces  are  planes  of  constant  Z',  with  the  bottom  face  at  Z'  =  Z  =  0. 
The  gage  positions  are  defined  by  a  letter-and-nuaber  pair;  the  letter  denotes 
the  face  and  the  number  denotes  the  gage  position  on  the  face. 

Table  B-l.  GAGE  POSITIONS,  COMPLETE  SUMMARY 


Gage 

Number 

Face 

X*  (cm) 
(Depth) 

Y'  (cm) 
(Width) 

Z»  (cm) 
(Height) 

Experiment  | 

E.l 

E.2 

FI 

Front 

0.00 

17.15 

5.69 

No 

No 

F2 

Front 

0.00 

14.68 

2.84 

Yes 

Yes* 

F3 

Front 

0.00 

12.24 

S.69 

No 

No 

F4 

Front 

0.00 

4.90 

8. S3 

Yes 

Yes 

B1 

Back 

11.79 

17. IS 

S.69 

No 

No 

B2 

Back 

11.79 

14.68 

2.84 

Yes 

No 

83 

Back 

11.79 

12.24 

S.69 

No 

No 

B4 

Back 

11.79 

4.90 

8. S3 

Yes 

Yes 

$1 

Leeward  Side 

S.89 

19. SB 

S.69 

Yes 

Yes 

S3 

Windward  Side 

S.89 

0.00 

5.69 

Yes 

Yes 

T1 

Top 

S.89 

14.68 

11.38 

Yes 

Yes 

T2 

Top 

S.89 

4,90 

11.38 

Yes 

Yes 

•High  frequency  response  ST-4- 11 6  gage 

All  gages,  except  for  the  g*ge  at  position  F2  for  shot  E.2  were  PcB 
Model  11SA24  gages,  with  a  frequency  response  of  SOG  kHx  and  a  sensitive 
element  diameter  of  0.S1  cm  (0.20  in).  The  gage  at  position  F2  for  shot  E.2 
was  a  Susquehanna  Instruments  Model  ST-4- 116  gage,  with  a  frequency  response 
of  2.0  HHz  and  i  sensitive  element  diameter  of  0.16  cm  (0.065  in). 


The  output  from  the  PcB  pressure  gages  was  recorded  on  magnetic  tape. 

The  tape  recorder  has  an  80  kHz  bandwidth,  which  provides  a  response  to  a 
square-wave  input  signal  to  within  4  percent  in  10  ns  and  1  percent  in  25  ps. 
Thus,  the  response  of  this  measuring  system  was  limited  by  the  tape  recorder’s 
characteristics  rather  than  by  the  gage  response.  The  digitizing  rate  for 
these  data  was  one  datum  every  2.5  ps. 

The  output  from  the  ST-4-116  gage  was  recorded  on  both  magnetic  tape 
and  an  oscilloscope.  The  response  of  the  ST-4-116/oscilloscope  recording  was 
limited  by  the  time  it  takes  the  shock  to  cross  the  gage.  The  ST-4-116/ 
oscilloscope  combination  was  used  to  obtain  data  as  close  in  time  as  possible 
to  the  passage  of  the  shock  front. 

The  digitized  data  were  analyzed  and  stored  on  the  BRL  CDC  173  for  later 
comparison  with  the  hydrocode  computations. 
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I.  COMPUTATIONS  H.l  AND  H.2 

The  finite  difference  grids  used  for  3-D  computations  H.l  and  H.2  (see 
Table  3)  are  identical.  Each  contains  92,512  flow  field  cells,  with  49  cells 
in  the  X  direction  (depth),  59  cells  in  the  Y  direction  (width),  and  32  cells 
in  the  Z  direction  (height).  The  target  is  modeled  using  6,912  nearly-cubical 
rigid  cells,  with  16  equal  X-direction  cells  (AX  =  0.7366  cm),  27  equal  Y- 
direction  cells  (AY  =  0.7253  cm),  and  16  equal  Z-direction  cells  (AZ  =  0.7112 

cm) .  The  rigid  cells  do  not  directly  enter  into  the  computation  other  than 

to  provide  a  non-responding,  perfect -reflection  structure  in  the  flow  field. 

To  conserve  storage  space  in  the  computer  and  yet  have  the  boundaries  of  the 
computational  grid  as  far  removed  as  practical  from  the  shelter,  cell  sizes 
moving  away  from  the  shelter  are  increased  using  geometric  progression  factors 
of  1.07  (or  less)  for  cell-to-cell  size  increases.  Cell  vertex  locations  are 
shown  in  Table  C-l. 

Figure  C-l  shows  a  top  view  of  the  initial  shock  location  in  the  hydro¬ 
code  grid.  The  shock  is  denoted  by  line  DE,  where  E  is  its  intersection  point 

with  the  X  axis  and  D  is  its  intersection  point  with  the  Y  axis.  The  axes* 

origin  is  denoted  by  point  P.  The  shock  velocity  vector  is  parallel  to  line 
PB,  which  is  perpendicular  to  the  shock  front  DE.  The  leading  corner  of  the 
shelter  is  denoted  by  point  A,  and  the  angle  of  obliquity,  a,  between  the 
shock  and  the  shelter  front  face  is  52.5  degrees.  Line  AC  is  parallel  to 
DE;  line  BC,  the  initial  distance  of  the  shock  from  the  shelter  corner  A, 
is  2.6551  cm. 


II.  COMPUTATION  H.3 

The  high-resolution  finite-difference  grid  used  for  the  3-D  computation 

H.3  is  a  modified  version  of  that  used  for  computations  H.l  and  H.2.  The  cell 

sizes  hav.*  been  reduced  by  a  factor  of  2  in  each  direction.  The  X  =  0.0  and 

Y  *  o.O  grid  boundaries  have  been  moved  closer  to  the  shelter  model.  The 

X  =  X  ,  Y  =  Y  and  Z  =  Z  boundaries  have  been  moved  in  past  the  respec- 
max  max  max 

tivo  ends  of  the  shelter  model,  so  that  not  all  of  the  shelter  is  modeled  in 
the  grid.  The  regions  of  primary  interest  on  the  shelter  were  those  near  the 
simulated  gage  positions  F2  and  S2.  The  grid  was  designed  so  that  enough 
computational  time  was  simulated  to  compute  tho  peak  overpressures  at  these 
points.  The  grid  reduction  was  necessary  to  keep  computing  time  within 
acceptable  levels.  (Halving  cell  sizes  in  a  3-D  grid  for  an  explicit-time- 
step  code  increases  computer  time  by  a  factor  of  24,  all  other  things  being 
equal . ) 

The  finite  difference  grid  contains  32,850  flow  field  cells,  with  45  cells 
in  the  X  direction,  73  colls  in  the  Y  direction,  and  10  cells  in  the  Z  direction. 
Tho  partial  shelter  model  is  built  using  12,190  nearly  cubical  rigid  cells, 
with  23  equal  X-direction  cells  (AX  =  0.3683  cm),  53  equal  Y-dircction  cells 
(AY  =  0.36265  cm),  and  10  equal  Z-diroction  colls  (AZ  =  0.3556  cm).  The  cell 
sizes  for  tho  grid  are  shown  in  Table  C-2. 
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TABLE  C-l.  CELL  DIMENSIONS  FOR  COMPUTATIONS  H.l  AND  H.2 


Index 

Cell  Dimensions 
(cm) 

Depth, 

I  Index 

Width, 

J  Index 

Height, 

K  Index 

1 

1.7521 

1.9634 

.7112* 

2 

1.6598 

1.8373 

.7112* 

3 

1.5723 

1.7192 

.7112* 

4 

1.4984 

1.6088 

.7112* 

5 

1.4109 

1.5055 

.7112* 

6 

1.3365 

1.4088 

.7112* 

7 

1.2660 

1.3183 

.7112* 

8 

1.1993 

1.2336 

.7112* 

9 

1.1361 

1.1544 

.7112* 

10 

1.0762 

1.0802 

.7112* 

11 

1.0194 

1.0108 

.7112* 

12 

.9657 

.9459 

.7112* 

13 

.9148 

.8851 

.7112* 

14 

.8666 

.8283 

.7112* 

15 

.8209 

.7751 

.7112* 

16 

.7776 

.7253  1 

j 

.7112* 

17 

,7366 

.7253*  | 

.7712 

18 

.7366* 

.7253*  ' 

» 

.7617  ’ 

19 

.7366* 

.7253*  ’ 

.8158  ; 

20 

.7366* 

.7253 

A 

.8737  ! 

(Continued) 

*  Rigid  cell  within  co&putetional  shelter  aodel 


50 


TABLE  C-l.  CELL  DIMENSIONS  FOR  COMPUTATIONS 
H.l  and  H.2  (Continued) 


Index 

Cell  Dimensions 
(cm) 

Depth, 

I  Index 

Width, 

J  Index 

Height, 

K  Index 

21 

.7366* 

.7253* 

.9357 

22 

.7366* 

.7253* 

1.0022 

23 

.7366* 

.7253* 

1.0733 

24 

.7366* 

.72S3* 

1.1495 

25 

.7366* 

.7253* 

1.2311 

26 

.7366* 

.7253* 

1,3186 

27 

.7366* 

.7253* 

1.4122 

28 

.7366* 

.7253* 

1.5124 

29 

.7366* 

.7253* 

1.6198 

30 

.7366* 

.7253* 

1.7348 

31 

.7366* 

.7253* 

1.8580 

32 

.7366* 

.7253* 

1.9899 

33 

.7366* 

.7253* 

34 

.7366 

.7253* 

35 

.78S8 

.7253* 

36 

.8383 

.7253* 

37 

i 

.8942 

,7253* 

38 

.9540 

.7253* 

39  • 

1.0177 

,7253* 

40 

1.0856 

.7253* 

(Continued) 

+Rigid  cell  within  coapuutional  skelter  model 
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TABLE  C-2.  CELL  DIMENSIONS  FOR  COMPUTATION  H.3 


Cell  Dimensions 
(cm) 

Index 

Depth, 

I  Index 

Width, 

J  Index 

Height, 

K  Index 

1 

.63300 

.65915 

.35560* 

2 

.63300 

.65915 

.35560* 

3 

.59965 

.61680 

.35560* 

4 

.59965 

.61680 

.35560* 

5 

.56805 

.57720 

.35560* 

6 

.56805 

.57720 

.35560* 

7 

.53810 

.54010 

.35560* 

8 

.53810 

.54010 

.35560* 

9 

.50970 

.50540 

.35560* 

10 

.50970 

.50540 

.35560* 

11 

.48285 

.47295 

12 

.48285 

.47295 

13 

.45740 

.44255 

14 

.4S740 

.44255 

IS 

.43330 

.41415 

16 

.43330 

.41415 

17 

.41045 

.38755 

18 

.41045 

,38755 

19 

.38880 

.36265 

20 

.38880 

.36265 

21 

.36830 

.36265* 

22 

.36830 

.36265* 

23  -  4S 

.36830* 

.36265* 

46  -  73 

.36265*  ! 

- 1 _ 

'  ‘Rigid  cell  within  computational  shelter  Model. 
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Because  the  partial  shelter  model  extends  from  the  bottom  to  the  top  of  the 
grid, .this  computation  does  not  model  rarefaction  waves  which  would  normally 
move  from  the  top  of  a  face  to  the  bottom.  However,  the  desired  simulated 
time  is  short  enough  that  these  waves  do  not  arrive  at  F2  and  S2  until  well 
after  the  peak  overpressure  occurs  at  those  points. 


III.  COMPUTATION  H.4 

The  finite  difference  grid  used  for  2-D  computation  H.4  is  a  2-D  slice 
of  the  grid  for  computations  H.l  and  H.2,  taken  at  the  ground  level  (K  =  1) 
plane.  The  columns  for  the  I  and  J  index  cell  sizes  in  Table  C-l  are  those 
for  the  computational  grid  for  H.4,  having  49  cells  in  the  X  direction,  and 
59  cells  in  the  Y  direction. 


.] 

i 

■  i 
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IV.  COMPUTATION  H.5 

The  finite  difference  grid  used  for  computation  H.5  has  twice  the  res¬ 
olution  (i.e.,  50  percent  smaller  computational  cell  sizes  in  each  direction) 
of  computation  H.4.  As  was  done  for  3-D  computation  H.2,  a  portion  of  the 
rigid  shelter  was  placed  in  the  upper  right  corner  of  the  grid  (the  high  I  and 
J  index  region) .  Enough  of  the  shelter  was  modeled  so  that  reliable  computa¬ 
tions  for  the  pressure  peaks  at  positions  F2  and  S2  were  possible.  The  compu¬ 
tational  cell  sizes  are  shown  in  Table  C-3.  The  cell  sizes  at  the  rigid  shelter 
are  AX  =  0.36630  cm  and  AY  »  0.36265  cm. 


V.  COMPUTATION  H.6 

The  finite  difference  grid  used  for  computation  H.6  has  twice  the 
resolution  of  computation  H.5,  and  four  times  that  of  H.4.  Its  configuration 
is  similar  to  that  for  H.5.  The  computational  cell  sizes  are  shown  in  Table 
C-4.  The  cell  sizes  at  the  rigid  shelter  are  AX  »  0.18415  cm  and  AY  »  0.181325  cm. 


VI.  COMPUTATION  H.7 

The  finite  difference  grid  for  computation  H.7  has  twice  the  resolution 
of  computation  H.6,  and  eight  times  that  of  H.4.  Its  configuration  is  similar 
to  those  of  H.S  and  H.6,  except  that  the  boundaries  have  been  moved  in  closer. 
The  computational  cell  sizes  are  shown  in  Table  C-5.  The  cell  sizes  at  the 
rigid  shelter  are  AX  »  0.092075  cm  and  AY  *>  0.0906625  cm. 
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TABLE  C-3.  CELL  DIMENSIONS  FOR  COMPUTATION  H.5 


Index 

Cell  Dimensions 
(cm) 

Depth, 

I  Index 

Width, 

J  Index 

1-6 

.77460 

.96400 

7 

.77460 

.$6315 

8 

.77460 

.94000 

9 

.77460 

.89000 

10 

.77460 

.85000 

11 

.77460 

.80900 

12 

.77460 

.76900 

13 

.77460 

.73200 

14 

.77520 

.69600 

IS 

.74600 

.66200 

16 

.70900 

.62900 

17 

.67500 

.59900 

18 

.64100 

.56900 

19 

.61000 

.54200 

20 

.58000 

.51500 

21 

.55100 

.48800 

22 

.52400 

.46600 

23 

.49800 

.44300 

24 

,47400 

,42100 

25 

.45100 

.40109  ; 

26 

.42900 

.38125  | 

r 

27 

.40700 

.56265  ‘ 

28 

.38700 

.36265  i 

29-30 

.36830 

.36265 

31  -  32 

.36830 

.36265* 

33-64 

.36830* 

.36265* 

65  -  64 

.36265* 

•Rigid  cell  vithin  coeputational  shelter  model, 
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TABLE  C-4.  CELL  DIMENSIONS  FOR  COMPUTATION  H.6 


Index 

— 

Cell  Dimensions 
(cm) 

Depth, 

I  Index 

Width, 

J  Index 

1  -  2 

.553125 

.633575 

3-4 

.553125 

.605 

5  i  6 

.575 

7-8 

.553125 

.545 

9-10 

.53 

.52 

11  -  12 

.505 

.495 

13  -  14 

.48 

.47 

IS  -  16 

.455 

.445 

17  -  18 

.435 

.425 

19-20 

.4125 

.4045 

21  -  22 

.3925 

.3845 

23  -  24 

.373 

.366 

25  -  26 

.3545 

.348 

27  -  28 

.3375 

.331 

29  -  30 

,3205 

.3145 

31  -  32 

.305 

.2995 

33  -  34 

.29 

.2845 

35-36 

.2755 

.271 

37  -  38 

.262 

,257S 

39  -  40 

.249 

.244 

41  -  42 

.237 

.233 

43  -  44 

.  2255 

,2215 

45  -  46 

.2145 

.2105 

47  -  48 

.2035 

.2005 

49  -  50 

.1935 

.190625 

51  -  58 

.18415 

.181325 

59  -  122 

.18415* 

.181325* 

123  -  166 

.181325* 

Rigid  cell  within  computational  shelter  model. 
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TABLE  C-5.  CELL  DIMENSIONS  FOR  COMPUTATION  H.7 


Index 


Cell  Dimensions 
(cm) 

Depth, 

I  Index 

Width 

J  Index 

.508375 

.5150125 

.4854 

.4908 

.4628 

.4671 

.4414 

.4451 

.4209 

.4242 

.4042 

.3828 

.3852 

.3650 

.3670 

.3481 

.3498 

.3319 

.3333 

.3165 

.3176 

.2019 

.5027 

.2879 

.2884 

.274S 

.2748 

.2618 

.2619 

.2496 

.2496 

.2381 

.2378 

.2270 

.2266 

.2165 

,2160 

1  .2064 

i 

.2058 

\  .1969 

.1961 

i  . 1877 

.1869 

.1790 

.1781 

.  1707 

.1697 

.1628 

.1617 

_ _ _ 

(Continued) 


Depth, 
I  Index 


1553 

1480 

1412 

1346 

1284 

1224 

1168 

1113 

1062 

1012 

0966 


Width, 

J  Index 

.1541 
.1468 
.1399 
.1333 
.1271 
.1211 
.1154 
.1099 
,1048 
.0998 
.0951 
.090662S 
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